Ultrafast spectral diffusion measurement on nitrogen vacancy centers in 
nanodiamonds using correlation interferometry 
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Spectral diffusion is the phenomenon of random jumps in the emission wavelength of narrow 
lines. This phenomenon is a major hurdle for applications of solid state quantum emitters like 
quantum dots, molecules or diamond defect centers in an integrated quantum optical technology. 
Here, we provide further insight into the underlying processes of spectral diffusion of the zero 
phonon line of single nitrogen vacancy centers in nanodiamonds by using a novel method based 
on photon correlation interferometry. The method works although the spectral diffusion rate is 
several orders of magnitude higher than the photon detection rate and thereby improves the time 
resolution of previous experiments with nanodiamonds by six orders of magnitude. We study the 
dependency of the spectral diffusion rate on the excitation power, temperature, and excitation 
wavelength under off-resonant excitation. Our results bring insight into the mechanism of spectral 
diffusion and suggest a strategy to increase the number of spectrally indistinguishable photons 
emitted by diamond nanocrystals. 



In the last decade the negatively charged nitrogen va- 
cancy center (NV) in diamond has emerged as a promis- 
ing resource for future quantum technologjE^. NVs, 
consisting of a nitrogen impurity atom with an adjacent 
vacancy, occur naturally in bulk diamond or diamond 
nanocrystals and are commonly used as single photon 
sources, operating even at room-temperature^. Diamond 
nanocrystals with a single NV are of particular inter- 
est, as they can be deterministically integrated into pho- 
tonic and plasmonic hybrid devices and circuits^S. The 
spectrum of the NV single photon emission consist of 
a broad phonon sideband around 700 nm and a narrow 
zero phonon line (ZPL) at about 638 nm. At cryogenic 
temperatures the ZPL of single NVs in bulk type Ha dia- 
mond has been shown to be almost Fourier-limited under 
off-resonant excitation, which allowed for the observation 
of two-photon interference | 12 * 13 1 , a crucial step towards 
applications in optical quantum information processing. 
However, for implanted NVs located close to the surface, 
several GHz wide spectral diffusion can be observed un- 
der off-resonant illumination!^ In nitrogen rich type lb 
nanodiamonds the energy levels fluctuate much stronger, 
leading to a broadened ZPL with a width of about 0.5 nm. 
It is widely assumed, that line broadening of emitters in 
condense d pha se is due to a fluctuating electrostatic en- 
vironmentl 15 ! 16 !. I n diam ond nanocrystals it is caused by 
ionized impuritie d 17 * 18 ! and charge traps. The spectral 
shift caused by single charges can be estimated using a 
sim ple toy model (sketched in Fig. la). Using data from 
RefsP^M] we found that the Stark-shift of the excited 
states A£ x / y — d ■ E in an NV center induced by a single 
elementary charge in a distance of about 10 nm can be as 
large as several hundred GHz. The time after which the 
probability that the electrostatic environment remains 
unchanged is reduced to 1/e can be defined as the spec- 
tral diffusion time T£>. Within tjj the system can be as- 
sumed to be free from spectral diffusion, and all photons 
emitted within tjj should be nearly indistinguishable. 
Enhancing the number of indistinguishable subsequent 



photons is a key goal. In this paper, we present mea- 
surements of td in NV centers in nanodiamonds under 
off-resonant excitation with a time-resolution exceeding 
previous experiments 21 with nanodiamonds by about six 
orders of magnitude. Our studies identify the excitation 
laser as the main source for enhanced spectral diffusion. 

The most straightforward method to measure spectral 
diffusion in single emitters is to take a time series of spec- 
tra and directly visualize the spectral wanderin g} 16 ! 22 ! 23 !. 
However, this method is only suitable for low spectral 
diffusion rates. For many emitters, like NV centers in 
nanodiamonds the rate of detected photons emitted from 
the ZPL transition is on the order a few kHZ, while spec- 
tral diffusion occurs much faster. 

In the literature two methods which in principle enable 
spectral diffusion measurements on timescales as fast as 
the emi tter's lifetime were suggested and demonstrated 
recently^4125l. The idea is to convert the frequency fluctu- 
ations due to spectral diffusion into intensity fluctuations. 
However, in this approach, technical limitations such as 




FIG. 1. (Color online) a) Toy model for spectral diffusion. 
Charge traps are ionized and lead to fluctuating electric fields 
at the NV center's location, b) Sketch of the experimen- 
tal setup. Milled type lb nanodiamonds are deposited on 
a solid immersion lens (SIL) and placed in a continuos flow 
He-cryostat. Single photon fluorescence is excited with a 
532 nm laser and collected through an objective lens with 
numerical aperture NA = 0.9. Using a spectral filter (BP) 
photons from the zero-phonon transition are filtered and sent 
through a folded Mach-Zehnder interferometer. 
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the required complex setup with several moving parts re- 
duced the time resolution in the first measurements to 
the order of 100 /ziP^l In a related approach Sallen et 
alP2l reached the subnanosecond regime, but for the price 
of a spectrometer-limited spectral resolution and an in- 
herently reduced detection efficiency. As pointed out in 
RefP3 this approach also requires the emission frequency 
to fluctuate around one fixed central wavelength, which 
is in general not the case for many quantum emitters. 
Here, we use an advanced and simplified interferomet- 
ric setup, combining the advantages of the previous ap- 
proaches. The timing resolution is limited only by the 
single photon counting instrumentation, which can be 
100 ps or below. The spectral resolution can almost be 
as high as for conventional Fourier spectroscopy, while 
the central wavelength of the emission does not need to 
be fixed. Furthermore, the photon detection efficiency 
is at least doubled compared to the scheme presented 
in RefPn Therefore, our method is applicable even for 
moderate integration times. 

We use a fixed Mach-Zehnder interferometer as dis- 
perse element, converting the spectral modulation of the 
ZPL into an intensity modulation which can be measured 
by correlating the photons from the two interferometer 
outputs. Figure lb shows a sketch of the setup. From 
the measured photon statistics, the timescale of the spec- 
tral diffusion is derived. In the following we calculate the 
cross-correlation between the outputs of the interferom- 
eter. We suppress the wavelength index A of the photon 
number operator h\{t) for simplicity. 
The time averaged photon detection rate in the left or 
right output port of the interferometer is < I{t)L/R >u 
with 



l(t) L / R = i] L/R h(t) m L / R (t), 



(1) 



where rj L / R is the overall quantum efficiency in the 
left and right exit of the interferometer, respectively. 
m L/n{t) is the interferometer introduced modulation 



m L/R(t) = lie sin[27rx/A(i)], 



(2) 



where c is the contrast of the interference fringes and 
x the path length difference of the interferometer arms. 
The plus and minus signs correspond to the left and right 
arm, respectively. The cross-correlation between the two 
arms ff^(r) = <^MW^);>* rea ds 



<lL(t)> t <I R (t)> t 



9lr(t)=9 {2) (t) < m L (t)m R (t + T) > u 



(3) 



where flf' 2 ^(r) = 



_ <.n(t)n{t+r).> t . g SCCO nd order au- 
tocorrelation function of the bare emitter and <: ... :>t 
denotes normal ordering before time averaging. 
We assume, that spectral diffusion occurs in form of 
jumps of the narrow emission line to random positions 
within a broad envelope^ and that individual jumps lead 
to a significant change of the interferometer transmis- 
sion. To calculate < mi,(t)m R (t + r) > f the character- 
istic function X(t,r), which is one if no spectral jump 



occurred within the time interval (i, t + t) and zero else 
is introduced and terms with and without spectral jump 
are sperated in the average. By using that A(t) = A(i+r) 
if no jump occurred and A(t + r) = A(i') ^ X(t) if a jump 
occurred, separating X{t,r) and identifying the proba- 
bility that the spectral position remains unchanged after 
the time r as p(r) =< X(t, r) > t we get: 

< m L (t)m R (t + r) > t = p(r) < m L (t)m R (t) > t 

+ [1 - P (t)\ ■ < m L (t)m R (t') > ttt , . (4) 



If the interferometer is adjusted in a way that several 
fringes are within the inhomogeneous width of the ZPL 
the sine-terms average out and < mz,(f) • m R (t') >t,t'= L 
while < rriL(t)m R (t) > t = 1 — c 2 /2. Therefore 



< m L (t)m R (t 



> t = 1-c 2 /2- P (t). 



(5) 



Remarkably, here also path length fluctuations on 
timescales slower than the SD average out and intcrfcro- 
metric stability is not required. Finally, the probability, 
that the spectral position of the ZPL remains unchanged 
after some time r is 



At) 
9 (2) {r) 



p(r) = 2/c 2 



1 - 



(6) 



Thus, it is sufficient to measure the autocorrelation 
<7^(t) of the bare emitter and the cross-correlation be- 

(2) 

tween the two outputs L/R of the interferometer g LR (r) 
to gain knowledge on the spectral dynamics of the emit- 
ter. Nanocrystals milled from high-quality type lb bulk 
diamond with a size of 30 nm to 100 nm were spin coated 
on a solid immersion lens (SIL) made of zirconium diox- 
ide^. The SIL is placed in a continuous flow He cryostat 
at about 5 K. Fluorescence is excited and collected in a 
confocal configuration (see Fig. lb) through a commer- 
cial NA 0.9 objective lens (Mitutoyo) placed inside the 
isolation vacuum of the cryostat. The NVs are excited 
by a green 532 nm laser with a power of several //W in 
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FIG. 2. (Color online) Spectra of the ZPL from a typical nan- 
odiamond at a temperature of 5 K with a removable band- 
pass filter centered at the ZPL (637 nm). The green curve 
(upshifted by 0.3 for clarity) is the NV center's fluorescence 
measured without the interferometer. The lower red curve is 
recorded after the interferometer. The inset shows the full NV 
spectrum with the yellow shaded area indicating the transmis- 
sion of the bandpass filter. 
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FIG. 3. (Color online) a) Measured second-order autocorre- 
lation function </ 2 '(t) (red line) and cross-correlation g^ R (r) 
(blue dots) of NV-1 at an excitation laser power of 14.1 /iW. 
The bin size is increased for longer delays r, to decrease the 
error. The error of </ 2 '(t) is in the order of the linewidth. 
The dip at zero time delay shows the single photon charac- 
ter of the emission, b) The probability p(r) that the ZPL 
remains constant on the order of 100 GHz within the time 
interval r calculated from the data shown in Fig. 3a accord- 
ing to Eq. (6). The bin size is increased for longer delays r, 
to decrease the error. The black line is an exponential fit to 
the data showing a spectral diffusion time td of (4.6±0.6) jis, 
while the dashed black line indicates the spectral diffusion 
time td (1/e time). 
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FIG. 4. (Color online) a) Dependency of the spectral diffusion 
rate on the excitation power, measured on NV-1. The blue 
curve is a linear fit to the data with the jump rate at zero 
excitation power being a free fit parameter. The red curve 
is a quadratic fit to the data. Obviously, the linear function 
fits better to the data, b) The number of collected photons 
from the ZPL transition per spectral jump as a function of 
the excitation photon energy, as measured on NV-2. The blue 
data points were measured at equal photon count rates and 
varying excitation powers (left axis). Point size corresponds 
to the error bar. The blue curve is a guide to the eyes. The 
red points show the spectral jump rate normalized to the ex- 
citation power, calculated from the same data set (right axis). 
It is clearly visible, that the spectral diffusion rate increases 
dramatically above the threshold of about 2.3 eV. 



front of the vacuum chamber. The single photons emitted 
by the NV are collected through the same objective lens 
and separated by a dichroic mirror from the excitation 
beam. An additional longpass filter blocks the excitation 
laser. In order to suppress the phonon side band, a re- 
movable bandpass filter (width 7 nm) centered at 637 nm 
is used. The filtered photons are directed into a Mach- 
Zchndcr interferometer and detected by two avalanche 
photo diodes (APD) in the two outputs and counted with 
a time-correlated single photon counting module (Pico- 
Quant). The arm length difference of the interferometer 
was adjusted to obtain 4 fringes per nm at 639 nm, the 
ZPL wavelength (see Fig. 2). To check the alignment 
of the interferometer one of the APDs was replaced by 
a 500 mm spectrograph (Acton SpectraPro 500i) with 
a cooled CCD detector and spectra were recorded (see 
Fig. 2). 

In order to derive the probability p(r) that the spectral 
position remains unchanged after the time r according 
to Eq. (6) first the bare second-order autocorrelation of 
the NV emission has to be measured. In order to do this, 
the whole fluorescence from the NV including phonon 
side-bands is sent into the interferometer. In this case 
the autocorrelation function g^> (r) measured by the two 
APDs is unaffected by the interferometer as proven by 
an independent measurement. The setup acts as a usual 
Hanbury Brown and Twiss setup. Next, the 637 nm 
bandpass filter is added in front of the interferometer 
and the cross-correlation g L ^(r) between the interferom- 
eter arms is measured with an typical integration time 
of 30 minutes. In this case the interferometer converts 
the spectral jumps of the narrow ZPL emission into in- 
tensity fluctuations. Care was taken not to change the 
laser power or adjustment within the total measurement 



time. 

The measured g^- and g^-functions of NV-1 taken at 
an excitation laser power of 14.1 fiW are shown in Fig. 3a. 
In order to derive the spectral diffusion probability p(r) 
the data is evaluated according to Eq. (6) , resulting in the 
data shown in Fig. 3b. Measurements taken with differ- 
ent interferometer adjustments revolving spectral jump 
widths from 260 GHz to 20 GHz showed, that the as- 
sumption of a narrow line jumping in a broad envelope is 
well justified (see online supplementary). Obviously, the 
integration time of 30 min was sufficiently long to resolve 
the spectral diffusion, but if higher timing resolution is 
needed the integration time must be increased. The de- 
crease of the probability p(r) fits well to a single exponen- 
tial decay. This supports the assumption that spectral 
diffusion is indeed caused by uncorrelated charge fluctu- 
ations. The spectral diffusion time tjj obtained from the 
fit is 4.6 ± 0.6 fj,s corresponding to a spectral diffusion 
rate jd = 1/t~d of about 220 kHz, while detection rate 
at the APDs of photons from the ZPL transition is only 
2.3 kHz. 

The achieved contrast c derived from the fit is 35%, which 
is lower than the measured interferometer contrast of 
90% due to fluorescence background and the two non- 
degenerate NV dipole transitions^. 
To get further insight into the origin of the spectral diffu- 
sion, we measured the dependency of the spectral diffu- 
sion rate 70 on the excitation power on NV-1. Therefore, 
the above described measurement was repeated for sev- 
eral laser powers between 3 /jW and 23 nW, well below 
the saturation power. The resulting data (Fig. 4a) show a 
clear linear dependency, a feature that was also observed 
on NV-3 where we repeated the measurement (see online 
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supplementary). Thereby we can rule out two photon 
processes like the photo induced charge conversion^! to 
be the origin of spectral diffusion. Remarkably for zero 
excitation power, the spectral diffusion rate reaches zero 
within the precision of the fit, although the intersection 
with the axis is a free fit parameter. From this we make 
two conclusions: First, the excitation with the green laser 
is the main cause of spectral diffusion. This is consistent 
with our simple charge trap model. The laser ionizes im- 
purities, providing floating charges which can be trapped 
in other charge traps. Second, the number of collected 
photons in the ZPL collected per spectral jump is con- 
stant for excitation far below saturation, i.e. independent 
of the excitation power. This number directly evaluates 
the quality of the single photon emission in terms of the 
possible number of subsequent indistinguishable photons 
available. This is a key figure of merit for future quan- 
tum optics experiments with NV-centers. Obviously, in 
the linear regime of excitation power (far below satu- 
ration) the quality of single photon emission cannot be 
improved by reducing the excitation laser power, as it is 
often done in experiments with self-assembled quantum 
dotsP. 

In order to derive a strategy to increase the number of 
collected photons per spectral jump we investigated its 
dependency on the temperature and excitation energy. 
Within the temperature range of 5 K to 20 K we did 
not observe any change in the spectral diffusion rate on 
NV-1. This is again consistent with the laser being the 
main cause of spectral diffusion by ionization of charge 
traps since their ionization energy exceeds k^T at this 
temperature range. 

In order to measure the influence of the excitation energy, 
we replaced the 532 nm cw laser with a pulsed super 
continuum source (NKT Photonics) with exchangeable 
10 nm broad bandpass filters. While keeping the count 
rate in the ZPL constant at (1.4 ± 0.2) kcts/s (well be- 
low saturation) we measured the spectral diffusion rate 
of NV-2 for several excitation wavelength from 510 nm 
to 580 nm, corresponding to photon energies between 
2.1 cV and 2.4 eV. We plotted the number of collected 
photons in the ZPL per spectral jump and the spectral 
jump rate normalized to the excitation power. The mea- 
surement in Fig. 4b clearly indicates, that the number of 



collected photons per spectral jump decreases and spec- 
tral diffusion rate increases, respectively, with increasing 
excitation energy. Remarkably, there is a pronounced 
threshold of about 2.3 eV. This gives evidence for the 
existence of deep charge traps with an ionization energy 
close to 2.3 eV. The remaining spectral diffusion for en- 
ergies below threshold might be attributed to substitu- 
tional nitrogen atoms in the diamond nanocr ystal s form- 
ing donor levels which are ionized at 1.7 eV * 17 * 18 *. More 
extensive studies to identify the trap states as well as to 
generally modify these states in order to reduce spectral 
diffusion and studies of higher quality diamond with nat- 
ural or implanted defects will be a subject of future work. 
Furthermore using two color excitation, the studies can 
be extended to energies below the ZPL transition to in- 
vestigate the influence of the nitrogen donors at 1.7 eV. 

In conclusion, we presented an interferometric method 
to measure the spectral diffusion of the ZPL of single ni- 
trogen vacancy centers in nanodiamonds. Our method 
works even if the spectral diffusion rate is orders of mag- 
nitudes higher than the photon detection rate and it is 
applicable for moderate integration times. Our system- 
atic studies show that in milled type lb nanodiamonds 
the jump rate is 1-2 orders of magnitude higher than 
the photon collection rate. Furthermore, the jump rate 
per photon is independent of the excitation power and 
NV temperature, but depends strongly on the excitation 
wavelength. Our results support the assumption that 
the main source of spectral diffusion stems from charge 
traps. A promising strategy to maximize the number 
of subsequent single photons from a nanodiamond is to 
use relatively high excitation powers, but clearly below 
the saturation level and to chose a proper low excitation 
wavelen gth. Actively enhancing the optical transition 
strengtrFED or collection rate is also required, possibly 
to enable active line stabilization schemes, as reported 
for NV centers in bulk diamond^. Finally, surface treat- 
ment accompanied by ultra-fast spectral diffusion studies 
should be extended to learn more about the dynamics of 
trap states on the diamond surface. 
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SUPPLEMENTARY INFORMATION 
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FIG. 5. (a): The spectral diffusion rate as a function of exci- 
tation power measured on NV-3. The blue curve is a linear fir 
to the data (green points). The red measurement points were 
not used for evaluation, as the count rate significantly changed 
during the measurement, (b): The single photon count rate 
as a function of excitation power, measured simultaneously 
with (a). The red points correspond to measurements, where 
the count rate was increased due to ablation of dirt. 



10 3 1? 10 5 10 6 10 3 It? 1? 10 6 

t / ns t / ns 

( c > 12 fringes/nm 





Exp(-T/5us) 

Measurement l 









io 3 it? io 5 io 6 



j I ns 

FIG. 6. Spectral diffusion time at different interferometer 
positions measured at NV-4. In (a) spectral jumps larger 
than 260 GHz are resolved. In (b) the resolution is on the 
order of 110 GHz, while in (c) jumps larger than 20 GHz are 
resolved. The spectral diffusion time for all measurements is 
the same, indicating, that spectral diffusion occurs in form of 
random jumps of the ZPL with a width below 20 GHz to a 
new position within the broadened ZPL. 



